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Using Wilson-Polchinski renormalization group equations, we give a simple new proof of 
decoupling in a 4 -type scalar field theory involving two real scalar fields (one is heavy with 
mass M and the other light). Then, to all orders in perturbation theory, it is shown that 
effects of virtual heavy particles up to the order 1/M 2N ° can be systematically incorporated 
into light-particle theory via effective local vertices of canonical dimension at most 4 + 2N . 
The couplings for vertices of dimension 4 + 2N are of order 1/M 2N and are systematically 
calculable. All this is achieved through intuitive dimensional arguments without resorting to 
complicated graphical arguments or convergence theorems. 
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1 INTRODUCTION 

The concept of effective field theories has played an important role in modern theoretical 
physics and it acquires its natural physical interpretation in the Wilson renormalization group 
formalism [3J. In the latter, one integrates out the high frequency modes scanned by a cutoff 
and then considers lowering the cutoff. This generates the RG flow, and the differential equation 
governing this flow is called the exact renormalization group equation or the flow equation. In 
general, as we scale down to a lower cutoff, the (Wilsonian) effective Lagrangian converges to- 
ward a finite-dimensional manifold parametrized by the "relevant" couplings which correspond 
to the renormalizable couplings in the conventional perturbation theory. This explains the phe- 
nomenological success of renormalizable quantum field theories in particle physics. Polchinski 
was able to give a simple proof of the perturbative renormalizability within this framework 
||, taking the 4 -theory as his model. Since the idea is very general, it has been extended to 
many different field-theory models, to study renormalization of composite operators and short- 



distance expansions, and also toward more practical applications. See for example p|— fllP . 

If we take seriously the above point of view on effective field theory and renormalization, 
it should be natural also to find, using this flow equation approach, an intuitive, yet rigorous, 



proof of the decoupling theorem [II] which states that, in a generic renormalizable quantum field 
theory with heavy particles of mass M, heavy-particle effects decouple from low-energy light- 
particle physics except for renormalization of couplings involving light fields and corrections of 
order 1/M 2 . In fact, as far as low-energy physics is concerned, heavy fields are not very different 
from massive regulator fields and hence those fields may be integrated out JT^| (together with 
high-frequency modes for light fields) as one lowers the cutoff below the heavy mass scale M. 
One then expects that the nonlocal irrelevant terms, generated in the effective Lagrangian by 
this integrating-out procedure, be suppressed as the cutoff is scaled further down, assuming that 
they are well-behaved. In this light, decoupling looks quite natural. Indeed, the decoupling 
theorem can be proved precisely along this line, and in this approach complicated graphical 
arguments or convergence theorems are entirely dispensed with. 

Since the resulting effective theory is renormalizable to the zeroth order in 1/M 2 , if one 
wishes to understand the low-energy manifestations of heavy particles, irrelevant (nonrenor- 
malizable) terms must be considered. This issue was investigated in Refs. |13|, fL4| and the 
result is that virtual heavy particle effects can be isolated via a set of effective local vertices 
with calculable couplings when combined with appropriate calculation rules to deal with irrele- 
vant vertices. This factorization theorem was proved to order 1/M 2 by adapting Zimmermann's 
algebraic identities |1| in the BPHZ formulation. The proof involves elaborate mathematical 



rearrangement and convergence theorems, and thus looks a bit artificial. [Also, as for the sub- 
graphs involving heavy and light particle legs simultaneously, the precise nature of factorization 
is not very clear in this approach.] In this paper, we present a rather elementary proof of this 
heavy- mass factorization to all orders in perturbation theory and to any given order in 1/M 2 , 
using flow equations. Actually the whole scheme allows a natural physical interpretation from 
the viewpoint of the RG flow. The calculation rules for irrelevant vertices are also systematically 
given. 
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There is, however, a subtle point. Our hope is to find a local effective theory with the 
UV cutoff A (^> M) which is accurate up to the order 1/M 2N ° (N is any fixed integer). In 
the original Wilson's view, the cutoff scale of the effective theory may well be identified with 
the heavy mass scale M, above which it is no more effective, and there is no need to worry 
about the presence of irrelevant terms in the Lagrangian in particular, since "natural scale" of 
those terms will be around M. Among them we may choose to keep explicitly some minimal 
number of irrelevant terms in our effective LagrangianF] for the accuracy of order 1/M 2N °. But, 
in the conventional discussion of quantum field theory, the UV cutoff A is supposed to go to 
infinity eventually. So, to connect it with Wilson's view, we may suppose scaling up the cutoff 
of the above Wilsonian effective Lagrangian from M to A . It will then generate infinitely 
many irrelevant terms which are unnaturally large. Also, during the scaling, all of them get 
mixed together and so we are forced to work with the Lagrangian consisting of infinitely many 
terms all the time. Any kind of truncation for the bare Lagrangian to some finite number of 
terms would yield divergences in physical quantities as A — > oo, because the unnaturally large 
coefficients would be amplified by some positive power of A /M as the cutoff is scaled down. 
This is nothing but the statement of nonrenormalizability in the language of RG flow. To avoid 
this problem we need to deal with irrelevant terms carefully, i.e., give suitable rules to obtain 
unambiguous finite results with only a finite number of terms included in the bare Lagrangian. 
It is achieved through the modification of the flow equation when irrelevant vertices are inserted, 
in the more-or-less same way as one treats composite operators and their normal products ]15| 



in the flow equation approach 0. Based on this procedure, we can connect the BPHZ approach 
with the present one. 

The plan of our paper is as follows. In Sec. 2, we present some basics on the flow equation, 
which are used throughout this paper. In Sec. 3, we define our model theory (the full theory) 
and choose a suitable form of cutoff functions which implement naturally the idea of effective 
field theory discussed above. We then obtain certain bounds on Green functions which are 
useful in the next section and prove the perturbative renormalizability. In Sec. 4, we prove the 
decoupling theorem and then proceed to show the factorization of virtual heavy-particle effects 
by deriving the local low-energy effective field theory which describes low-energy light-particle 
physics to any desired order in powers of 1/M 2 . We conclude in Sec. 5. In the Appendix, we 
briefly describe the renormalization of Green functions with single or multiple insertion of local 
vertices and Zimmermann's normal product which are needed in Sec. 4. 



2 THE FLOW EQUATION 

In this section we present some features of the flow equation in general setting to fix 
notations and to facilitate our later analysis. For simplicity we consider a theory of a single 

1 For example, the continuum version of Symanzik's improved action [ft6| in lattice theory have been discussed 
in Jl7| by adding suitable irrelevant terms in Lagrangian. 
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scalar field in four Euclidean dimensions with a cutoff Ao- The bare action is written as 

S^] = i(0,P^V) + i>[0], (2-1) 

where P A ° is the free-particle cutoff propagator, (/, g) is defined by the momentum-space 
integral of / and g, (f,g) = / j^zf{p)g{—p) , and L Ao [<f)} represents the interaction part of 
the bare action. We will also allow the insertions of some additional local vertices or certain 
composite operators (of unspecified physical origin). To account for this, we define Ltat[0] as a 
formal power series in a = («i, . . . , «&), which has L A °[0] as its zeroth order term, viz., 

L tt W\ = E lA °'' N W > ( LA °M = LAo; °M) ( 2 - 2 ) 

|JV|>0 ' 

where N is a multiindex N = (N u ...,N k ), \N\ = Ei=i N h N\ = iVi!---iV fc !, and a N = 
ai 1 ■ ■ ■ a^ k . L Ao;iV 's may be regarded as additional local vertices appended to the original 
Lagrangian L Ao or as composite operators in which one is interested. Also, we will denote 




|AT|>1 



The generating functional, with the insertion of the operator e L c°M, is 

Z[J] = J D0 e -^*' pA °" 1 ^- i totM+<^). (2.4) 

Following Wilson and Polchinski [fll 0], one may integrate out the high- momentum components 
of 4> and reduce the cutoff A to a lower scale A. The result is []5| 

Z[J\ = J p0 e -i<^-V>-iL[^° J +*]+< J .*>+3W^J) , (2.5) 

where P A (P^ ) is the low-(high-)frequency part of the propag atof| with P Ao = P A + Pf , and 

L A ot is the generating functional of amputated connected Green functions (with e L c^ inserted) 
with both UV cutoff A and IR cutoff A. Physical Green functions are obtained from L A ot in 
the limit A — ► Jl, ||, ||. On the other hand, if J(p) = for p > A so that J couples to the 
low-frequency modes only, all J's drop out from Eq. ( |2.5| ) except for (J, <p) since P^° has only 
high-frequency modes, and L A ot coincides with the Wilsonian effective action [|IJ . L A ot obeys the 
exact renormalization group equation or the flow equation [[|, 

2 There will be infinite ways in doing the separation; choosing one way of separation may be considered as 
choosing a "renormalization scheme" in the Wilson renormalization group approach, and physical quantities 
are independent of such choices. While it is not necessary to specify a particular scheme at this stage, we will 
choose a specific scheme later to facilitate our discussions. 



2 t A 



PL 



tot 



SL 



tot 



_5(t){j))5(j)(-p) 6<f>(p) 6cf>(-p) 



(2.6) 



fielH-den. nart 
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To fix appropriate boundary conditions for L A ot , renormalization conditions in the conventional 
renormalization procedure may be used at A = where physical Green functions are generated. 
Expanding L A ot in powers of a, we now define L A;Ar, s as the expansion coefficients, i.e., 



L 



totL 



a N 

\N\>0 



A\N r 



[0], (L Afi {<P]=L 



(2.7) 



Then, L A ' ,N [<j)} is the generating functional of amputated connected Green functions with an 
insertion of an operator On which is identified as the N-th order coefficient of e L c ^, i.e., 



v-^ OL N ^ L A 
\N\>1 ' 



ex P (E ^ A ° ;7V ) 

|7V|>1 iV • 



N 



(2- 



Explictly, On can be calculated as 



O 



N 



e ni n 



1 /L Ao;/ \ n/ 



|J|>1 



nj! 



(2.9) 



where J = (Ji, . . . , I k ) and J! = I x \ ■ ■ ■ I k \. For example, if k = 1, Oi = L Ao;1 , 2 = L Ao;2 + 
(L Ao;1 ) 2 , 3 = L Ao;3 + 3L Ao;2 L Ao;1 + (L Ao;1 ) 3 and so on. At the zeroth order of a, we obtain the 
flow equation for the effective Lagrangian L A from Eq. (|2.6| ), which assumes an identical form 
as Eq. (|2.6| ) other than the replacement L A ot — ■> L A . At order a N , on the other hand, we have 



d\ L A ' N 



d 4 p 



d A P\ P ) 



2 T A:N 



6 2 L 



5L A 5L A,N 



8<p(p)6<f>(-p) S(p(p) 6<p(-p) 



E 

0<KN 



V I J 8<j>{p) 8<t>{- P ) 



(2.10) 



field-dep. part 



where 



and I < N means Jj < A^j (for every 



k) while 



(7) -(£)-(£. 

I < N represents I < N and |/| < |AT|. If k = \N\ = 1, i.e., for a single insertion of operator 
through L Ao;1 , the last term vanishes and so we have a linear and homogeneous equation in L A;1 , 
which has been used in discussing composite operators @, || . Equation ( 2.10|) with k = 2 and 
N = (1,1) has also appeared in the literature || in deriving the short- distance expansion of two 
composite operators. Note that, if \N\ > 1, Eq. ( |2.10|) contains inhomogeneous terms. However, 
in spite of the inhomogeneous piece, it is still of first order in L A ' N while the homogeneous part 
remains exactly the same for all N. Therefore the general solution will be a sum of a particular 
solution and a solution to the homogeneous equation (which is just the 1 7V| = 1 equation for 
Green functions with a single insertion of an operator). This property will be used in the 
Appendix where Zimmermann's normal product and multiple insertion of local vertices are 
discussed in this formalism. Equation ( |2.1(J| ) is the starting point of our subsequent analyses. 
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3 STRUCTURES OF THE FULL THEORY 



3.1 Defining the full theory 

As our full theory we consider a scalar theory ((p-ip theory) which involves two real scalar 
fields and ip, defined with an ultraviolet cutoff A . The mass M of ip field is assumed to 
be much larger than the mass m of field but of course much smaller than A . In the next 
section, we will discuss the effect of the heavy field ip at low energy. The bare action for the 
system reads 

S (/)Ao 4<^°^) ' 2 

where P^° and P M ° are respective free-particle cutoff propagators. We assume that particles 
interact via quartic couplings, observing Z 2 x Z 2 symmetry. Explicitly, we may write 



(3.1) 



^(/)Ao 



/ 



f 

PA ,2 , P5 I o , \ 2 , P% 1 4 , Pi 1 2 1 2 



(3.2) 



The superscript / is used for quantities of the full theory. As is conventional, the bare couplings 
pi (a — 1, 2, • • • , 7) may be written as 



p[ 

A 



Z^m\ - m 2 , 
Z^M 2 - M 2 , 



p{ 
Pi 



Za — 1, 



z^ — 1, 



pi 
pi 



(3.3) 



pi — ZaZ^qI , 

where (mj^M, 2 ), (Za, Z^) and (g®, g®, g®) are bare masses, wave-function renormalizations and 
bare coupling constants, respectively. Propagators P^° and P M ° are defined by 



pA = Rm{Kp) 



p2 _|_ m 2 



pA 



R M (A ,p) 
p 2 + M 2 ' 



(3.4) 



where R m (A, p) and R M (A,p) are smoothed variants of the sharp cutoff function 6( A— p). As in- 
dicated in Sec. 2, choosing R m and Rm corresponds to choosing a "renormalization scheme" and 
physical quantities are not affected. Here we choose the following "mass-dependent scheme" : 



Rm(Kp) 

Rm(Kp) 



1 - if ((1 + A/m) 2 
1 - if (A 2 /M 



if 



k £ 



A 2 



A 2 



(3.5) 



where if is a C^-function on [0, oo) such that if (a) — 1, < a < 1 and if (a) = 0, a > 4 with 
if (a) — > 1(0) exponentially as a — > 1(4). Thus, both R m (A,p) and Rm(A,p) are C°°-functions 
on ([0, oo) x R 4 ) and i? m (i?M) approaches exponentially as A — > (as A — > M). If A > m for 
i? m or A > 2M for R M , they simply reduce to if (p 2 / A 2 ). A notable property of Rm is that if 
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A < M, we have Rm = identically. This choice of Rm is very natural for our purpose, because 
it implies that at the scale A = M we have all modes of the heavy field i\) integrated out and 
there remains only the light field below the scale M; it explicitly implements Wilson's point 
of view on effective field theory. Therefore, if A < M, the flow equation will look like that 
for the theory having the light field only. One may regard this property specific to our choice 
( p.5| ) as the analogy of the so-called "manifest decoupling" in the conventional approach |L8 . 
On the contrary, if one chose a "mass-independent scheme", i.e. if R m and Rm were chosen 
independently of their masses m and M, a substantial part of heavy particle modes would not be 
integrated out even below the scale M, thus making subsequent discussions rather complicated. 
Now, we define L^ A [(p,ip] following the general line discussed in Sec. 2. It will then satisfy 



dhL ~ 2 J (2tt) 4 \ dAjP ' 



(p) 



SL^ A 5L^ A 



6(/>(p)5(f>(-p) 5</>(p) 5(j)(-p) 



+ (m -> M, 



We expand in powers of and ip in momentum space, utilizing Z 2 X Z 2 symmetry, 

r 2ni rl 4 r>- 2n2 ~ l d A ri- 

£ (/)A [^]=£ £// II 7^ II 11 



field-dep. part 

(3.6) 



|r|>Q |n|>l 



A (2vr) 4 (2tt) 

x £>r£n (Plj ■ ■ ■ :P2ni'iPl: ■ ■ • >P2n 2 -l) 



0(Pl) • • • 0(P2m)^(pi) • • • ^(P2na-l) 



(3.7) 



where p' 2ri2 = —Y,f=S.Pj ~ X^S -1 Pj- Some explanations on our notations are in order: g = 
(91,92,93) are perturbation-expansion parameters which may be identified as renormalized 
coupling constants (see below), r and n represent (ri, ^,7*3) and (711,712), respectively, with 
\n\ = Hi + n>2, \r\ = r\ + r 2 + r% and finally, g r = g^g^g^- Then the vertex function £f2n nas 
Bose symmetry and, in addition, will possess following properties: (i) as we noted in Sec. 2, 



C^ln ° can be identified with amputated connected Green functions G^n °f the theory, i.e., 



(Pi) • • • )P2ni)Plj • • • ;P2n 2 -l) ~~ ^V^n(Pl> ■ ■ ■ iP2n 1 ]Pii ■ ■ • ;P2n 2 -l) ' 



(3i 



(ii) £,^2n = if |n| > |r| + 1 since just the connected diagrams contribute to £^2n as we have 



seen in Sec. 2; (iii) € C°°([0, A ] X [R 4(n 1 - ) ) and is invariant under the 0(4) rotation group. 

For later use, we also denote vertex functions obtained by summing over r 2 and r 3 as C, i.e., 



Mm 

*~Ti,2n 



V n r2 n r3 r (/)A 

i/2 i/3 i '(ri,r2,r 3 ),2n ' 



(3.9) 



We insert the expansion ( 377 ) into the flow equation 

d A^rln (Pi, P2nx 5 Pi, • • • , P 2 , 



2 n 2 -l/ 



2m + 2 

2 



(27r) 4 ,2 ^ ni+1,n ^ 



and arrive at the equation 



,P2m,g, -g;p'i,---,p 2 „ 2 -i; 
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+2 \ hl[d A P^P) [C^ hM {p li ... iPaiL ^P ] ^... f ^ h _ 1 ) 

<P2ti 2 -1j 
,P2ni,Pl> ■ ■ ■ ?P2n 2 ' <l) 

«(/)A /-„ „ .„/ „/ 



; 2+ ! 2 = n 2 

_ /2n 2 + 2 
I 2 

+2 E 



(2vr) 4 



symm 



z: (/)A r», 



r' +r" — r 



( 2 +Z^ = n 2 + l 



X A-",2(Zi,^)(P2Zi+l, • • • ,P2 ni ',P2l 



2' - - - ' 



P2r, 



symm 



(3.10) 



where P = - J^fli 1 Pj - E|=i Pj, P' = - Ej=i Pi - E^f 1 Pj an d [" " ijmm implies symmetriza- 
tion with respect to momenta Pi, ■ ■ ■ ,P2m (and momenta p[, . . . ,p' 2n2 )- Denoting 0(^)-legs by 
thin(thick)-lines and d\P A by a straight line, we can represent this equation by the diagram 
shown in Fig. 1. Here some remarks on Eq. ( p.lOj ) are in order. First, since we have C^lt = 
(see below) the restriction r' + r" = r also implies that r',r" < r (i.e., \r'\, \r"\ < \r\ and 
r 'i-, r i < r i)- Also, in relation with the effective theory of the next section, it should be noted 
that, for Cfflr ni o) (i.e., no external heavy particles), the last term is identically zero and ip 
enters the flow equation only through the A-differentiated propagator in the third term. 

If suitable boundary conditions are given, the flow equation (|3.10|) will produce the ampu- 
tated connected Green functions of the theory, G r2n 



r (f)A=o 

*-'r,2n 



. From Eq. ( |3.2|) , we have 



A = A, 



o • 



aw />(/) A 
u p J ~'r,2n 



0. 



d™£ 



(/)Ao 
r,2n 



(0) ~ bare couplings in Eq. (3.3), 
(to be determined) 



2\n\ + \w\ > 4. 
2\n\ + \w\ < 4 . 



(3.11) 



[Wjl, 



where w is the multiindex {w} = {wi, . . . , w 2ni , w[, . . . , w 2 n 2 -i}> w j 
d™ 1 ■ ■ ■ Cn" 1 ^ • • • d^h , = nl=i d w ^/{dp^) w ^. We will also write \w\ 



and 



'Pin, 1 "Pi ~ "Fi " ' \--f'F/ ' "- i~i — ^,/i\ wl C 1 

At A = 0, on the other hand, we impose following renormalization conditions on relevant 
terms0: 

*(/)A=0^_ n g d n „£ im=0 < 



£ 



£ 



r,(2,0) 
(/)A=0 / 



(0) = o, 



r,(0,2) (Pi 
£(/)A=0 



0. 



J Pn u Pv'~'r,(2,0) 
u Pp U Pv ^,(0,2) \Pl, 



'(o) = o, 



£ 



(/)A=0 



0. 



r,(4,0) 
.(/)A=0 



(0) 



4{Or,(l,0,0) > 



A", (0,4) (P2)P3'P4) 



41^,(0,1,0), (3.12) 



,,,2.2) ( 'P5) = 1^,(0,0,1)- 



Here, normalization momenta p\ (i = 1, . . . , 5) are chosen to be constants of magnitude M; 
i.e., we have chosen the renormalization points for the light-particle Green functions at zero 

3 With these conditions the expansion parameters (51,32,33) are now identified with renormalized coupling 
constants. As far as our proof goes, details of the renormalization conditions are not important as long as 
one choose the renormalization points for the light particle at momenta much smaller than M fl8|| . But, for 
practical applications, it would be necessary to choose the renormalization conditions for the heavy particles at 
momentum values of order the heavy particle mass pjj] . 
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momentum, and those for the heavy particles at momenta of magnitude M. It is then easy to 
show that p{ are uniquely fixed within perturbation theory. (See, for example, Ref. 0.) 

3.2 Bounds on Green functions 



Now the theory being completely specified, one may use the flow equation to demonstrate 
the perturbative renormalizability |2], |5|, which amounts to boundedness and convergence of 
the vertex functions £^2n as ^o ~~ * 00 • To reach our goal in the next section, however, we need 
slightly stronger bounds which take particle-mass dependence into account explicitly 

First, we introduce a set of norms on vertex functions || ||( a ,6,b'), for positive real a, b, and b', 



f)z r (f)A 

U i ~r,2(n 1 ,n2) 



(0,6,60 



(3.13) 



|to|=« 



For vertices with no heavy-particle leg, we will often omit b', 
Now it is not difficult to show that, with the choice ( |3.5| ), we have the following uniform bounds: 



xzAf)A 



(o,6) 



az r{f)A 



(a,b,b') 



d z p d A R m (A,p) 



(p 2 + m 



2\n 



< 



(A + m) 2n+z+1 



d z p d A R M (A,p) 



( P 2 + M 2 y 



< 



^271+2 + 1 



6 (A - M) 



(3.14) 



where < A < A and c, d are constants which do not depend on A, A , M, or m (but may 
depend on z, n and k). Acting on Eq. ( 3.1 0|) and estimating the resulting expression with 
the help of the bounds ( |3.14j ), we find that, for a fixed constant /i of order m, 



a az Af) A 



(2A,/i,M) 



(3.15) 



< const | (A + m) 



U ^,2(711+1,712) 



(2A,/i,M) 



+ 0(A- M) 



azAf)A 

u / -7-,2(ni,n 2 +l) 



^ r / + ,.// =r r / +r // =r / (A + 171 



I 3+21 



p,z 2r U)A 



(2A,n,M) 



(2A,n,M) 



az a Af)A 

° L 'r",2{l' 1 ,l' 2 ) 



{2A,n,M) 



!+!'=(ni+l,n 2 ) 
z l+ z 2+ z 3= z 



i+i'=(™l>«2+ 1 ) 
z l+ z 2+ z 3= z 



where "const" stands for some finite number independent of M, A and Ao (but may depend on 
m and \i through (/i/m)). Equation ( |3.15| ) holds for all A G [0, Ao], and it allows us to present 
our subsequent results and proofs in simpler forms than those appeared in earlier works ||, [5|. 
Now we are ready to prove the boundeness of vertex functions: 



Theorem 1 For vertices with no external heavy-particle leg, 



azAf)A 

u ^,2(711,0) 



< (A + m) 4 -' ni ^Plog 

(2A,n,M) — v ; ta 



A + m 



m 



< A < A . 



while, for irrelevant vertices (2\n\ + z > A) with external heavy-particle legs, 



(f)A 
r,2n 



(2A,fi,M) 



< (A + m) 4 - 2|n| - z Plog 



A + M 



< A < A 



(3.16) 



(3.17) 
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and, finally, for relevant vertices (2\n\ + 4 < 4) with external heavy-particle legs, 



(/)A 
r,2n 



(2A,n,M) 



< (A + M) 4 - 2|n| - 2 Plog 



A + M 



m 



< A < A , 



(3.18) 



where Plog(-) is some polynomial in log(-) whose coefficients are independent of M, A, and A . 

Before proving Theorem |l], we here note that the above bounds are not unexpected — they 
just show the right scaling behaviors. That is, the effect of imposing renormalization conditions 
at momenta of order M for heavy-particle legs shows with appropriate powers of M (up to 
logarithmic corrections) in the flow equation of vertices with external heavy-particle legs. On 
the other hand, in the light-particle sector, the bounds have the same form as those of a single 
scalar theory [0, ||: the large M corrections do not show up because m is forced to be small by 
hand. Of course, this should be the case for decoupling to occur in the first place. 

Proof of Theorem We proceed along the induction scheme set up in Ref. [0]. Equations 
( p.l6| )-( p.l8| ) are trivially true for r = or for (r, n) with |r| > and |n| > |r| + 1, because 
Ju/Jn = f° r those ( r i n)-values as we have seen before. Now suppose it holds for (r, n) satisfying 
the condition |r| < s — 1, or |r| = s and \n\ > k + 1. Then the induction step consists in proving 
Eqs. ( |3.16| )-( |B.18D for any fixed (r, n) with |r| = s and |n| = k. We begin with the irrelevant 
vertices, i.e., with the case 2k + z > 4. For these we use the boundary condition ( |3.11|) to write 



azAf)A 
u 1 ~'r,2n 



< 



Mo 



(2A,m,M) - J A 



(3.19) 



The right hand side is now estimated with the help of Eq. ( 3.15 ). If A > M, according to the 
induction hypothesis it is easy to find that, for irrelevant vertices with no heavy-particle 



° ^r,2(ni,0) 



(2A,n,M) 



< 



/■Ao 

/ d\(X + m 

J A 



< (A + m 



A-2k- 




(3.20) 



A similar estimate yields the bound fl3.17| ) for irrelevant vertices with heavy-particle legs. For 



A < M, the bounds on d™ C[^ 2 ^ with z — 0, 1 in Eq. ( |3.18| ) look problematic because they are 
not suppressed by appropriate powers of A, i.e., they might result in undesirable extra powers of 
M when one estimates the right hand side of Eq. ( |3.15| ) with the induction hypothesis. They, 
however, cause no real problem since they appear in the right hand side of Eq. ( |3.15| ) only 
with 6 (A — M). There are other dangerous terms when one estimates vertices with n 2 = for 
A < M, namely, the log M/m corrections in Eqs. (|3.17| ) and (|3.18j ); they are again multiplied 
by 8 (A — M) and do not contribute at all if A < M. Therefore, we conclude that the induction 
argument used in A > M case goes through with no obstuction for irrelevant vertices. 



*Each time the symbol Plog appears, it stands in general for a new polynomial of log. 
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For the case of relevant vertices (2k + z < 4), we begin with 4^(04)' Recalling the boundary 



conditions at A = 0, we integrate the flow equation (|3.10| ) from to A (at p\ = p'^ to obtain 

r-A 



l4$>1)l < l4^VJI 



">,(0,4) 



dX 



r r (f)x 



(2A,fi,M) 



< Plot 



A + M 



m 



(3.21J 



where we have used the induction hypothesis in the second line. Now, using the Taylor's formula 



4^)(K) = 4?m)(p1) + E(ri-pJ) f dtd^c^wi^ 

i=i j ° 



(3.22) 



and then putting A = M,we get a bound for 4^(04) as being of order Plog(M/m) at most. 

Then the integration of ^4^4) fr° m A to M gives Eq. ( |3.18| ). The same argument works for 
other relevant vertices with 2k + z = 4. Next, if n = (0, 1) and z = 1, we may use the property 
that dpC^2n (0) = f° r °dd dimensional vertices (due to Z 2 x Z 2 symmetry), and therefore 



a>L >,(0,2) 



(2A,/Lt,M) 



< (A + M)Plog 



A + M 



ni 



(3.23) 



The same argument shows that 



BC U)K 

Ui ~r,{2,0) 



< (A + m)Plog fA±™y it i s then straight- 

forward to obtain the bounds for \n\ = 1 and z = using Taylor expansion around the 
renormalization point, as we did for 2k + z = 4 case. I 

In particular, at A = 0, Theorem 1 tells us that amputated connected Green functions and 
their derivatives of the theory are bounded, for \pi\ < fi and \p' { \ < M, by 



const • m 4 2ni |w| , n 2 = 

I^g^Cpi, • • ■ , iwi, • • • ,pU-i)I < { - 4 - 2|n|HHpi °g (f ) > ^ ^ o. 2 I-I + H > 4 

M 4-2|n|-M plog (M) ; n2 ^ o, 2|n| + M < 4. 

(3.24) 

Now we show that Green functions converges to a finite limit as the UV cutoff A goes to 
infinity. It will also complete the proof of the perturbative renormalizability of the <p-ip theory. 



Theorem 2 For vertices with no external heavy-particle leg, 

'A + m\ 3 /K . sa _,„,__ Mo 



< 



(A + m) 3 - 2ni -*Plog (—) , < A < A , (3.25) 
V m J 



(2A,n,M) ~ V A 

while, for irrelevant vertices (2\n\ + z > 4) with external heavy-particle legs, 



<JA <J >~-r,2n 
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and, finally, for relevant vertices (2\n\ + 4 < A) with external heavy-particle legs, 

IK^IU,^ < (^) 3 (A + W)«W-Ptog (£) , < A < A . (3.27) 

Proof. The method is very similar to that employed in Theorem 1 and we only sketch the 
procedure. For irrelevant vertices, we use the identity 

dh dp(C^l n ° — C^n ) — dp^A^-l^n U=A + J d\dxdA dpC^n, (3.28) 
while, for relevant vertices, we use the relation 

/•A 

d Ao dp(C r2n - C r ,2n~ ) = J q d\d x d Ao dp C r<2n . (3.29) 

Then, differentiating Eq. (|3.10| ) with respect to Ao, we obtain the bounds for the right hand 
sides, as in Eq. fl3.15| ). Finally, applying the induction scheme as described in Theorem 1 and 
using the results of Theorem 1, it is straightforward to prove the theorem. I 



4 DECOUPLING, FACTORIZATION AND EFFECTIVE 
LAGRANGIANS 

We are now ready to discuss the low-energy physics of the (fi-ip theory. In this section we 
will prove the decoupling theorem and the factorization of virtual heavy-particle effects to any 
desired order in 1/M 2 on the basis of the flow equation approach. We start with proving the 
decoupling theorem which corresponds to the zeroth order result of factorization. 



4.1 Low-energy effective field theory and the decoupling theorem 

Let us consider a theory of light particles (the theory) of mass m with the interaction 
part L A ° of the bare action given by 



-Ao 



(Tic 



(4.1) 



As before, L A [(j)} is defined to satisfy the flow equation ( |2.6| ) and can be expanded as 



oo oo . 2n-l 74 . 



r=0 n=l 



,P2n-l; 



(4.2) 



The vertex function C A 2n has properties similar to those for C^ln'- £-r2n has Bose symmetry, 
C A 2n = if n > r + 1, L A 2n is C°° and Euclidean-invariant, and finally C A 2 ^ is equal to the 
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amputated connected Green function of the theory, i.e., G c r2n . It satisfies the flow equation 



dA^ 2n (p u ... J p 2 n-i) 



/2n + 2 
V 2 



d A q 



(2vr 



(4.3) 



+2 ]T ll ' d ^(P) Uv lJB (Pl, • • " ,P2l-l)C? ll>2l ,(-P,P2l, • • • ,P2n-l] 



r -f-r =r 
i + !'=n + l 



symm 



where P = — J2j=i Pj- The boundary conditions which C A 2n 's obey are determined as follows. 
First, at A = Ao, irrelevant vertices vanish, i.e., 

<9;<° 2 „ = 0, 2n+M>4. (4.4) 
At A = 0, we impose the same renormalization conditions ( |3.12| ) as those in the full theory, 

<F°(0) = d Pll d p „£^°(0) = 0, <I°(0) = ±6 rl , (4.5) 

so that the <fi theory may not be distinguishable from the (f>-ip theory at low energies. 

From Theorems [l] and |2| we immediately have the bounds on the vertex functions C A 2n 
(which amounts to proving the perturbative renormalizabililty of the <fi theory), viz., 

'A + m s 



(2A lA i) 



< (A + m 



A-2n-z 



Plog 



^ oz nA 
JA u t~ r2n 



(2A,M) 



< 



A + mV 



m 

\3—2n—z 



A r 



Plog 



< A < A 
A 



(4.6) 

0<A<A . (4.7) 



-J (A + m) c 

Now we are in a position to prove the decoupling theorem: 
Theorem 3 The difference of the vertex functions £^2(n0) ^ n theory and C A i 2n in the 



theory, where r = (7*1, r 2 , r 3 ) , satisfies the bound 

'A + m x 2 



° \'~'r,2{n,Q) z -ri,2n°r,(ri,0,0) ) 



(2A.A0 



< I 



M 

^A—2n—z 



(A + m) 4 - 2 "- 2 Plog 
Hog (A). 



< A < M 



M < A < A , 
(4.8) 

Proof. If A G [M, Ao], the bound is immediate from Eqs. ( |3.16| ) and ( [4.6| ), and hence it suffices 
to consider < A < M only. Let us define V A 2n as the difference of vertex functions, 



V A = C 



(/)A 



-a 



(4.9) 



-r,2(n,0) / ~'r 1 ,2n O r,{r 1 ,0fl) ■ 

Subtracting the flow equation ( |3.10| ) for Cj 2 t fr° m that for C A 2n (given in ( |4.3| )) and making 
estimates using Eq. (|3.14| ), one easily finds that, for < A < M, 



d A d z V^ 2n 



(2A, M ) 



< const < (A + m) 



9 Z ^r,2(n+l) 



E 



r +r — r 
i + ;'=n + l 
z l + z 2 + z 3= z 



(A + m) 3+2 



r',2/ 



(2A, M ) 



(2A,/n) 
f) Z3 T) A 

u u r t, 2V 



(4.10) 



(2A )A 



+ 



f)Z 3 pA 



(2A,„) 5r "'«'A0) 
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where it should be noted that terms involving the heavy-particle propagator all vanish. Bound- 
ary conditions for V^ 2n are simply 



f) w 7~) Ao - 

u p u r,2n - 

P r A 2 =°(0) 



0, for 2n + M > 4 



P A 4 =0 (0) = 0. 



(4.113 



Since V^ 2n 



0ifr = 0orn>|r| + l with |r| > 0, the theorem holds trivially for these 



values of r, n. So we will now prove the validity of Eq. 



assuming that it is true for r and 
= s. As usual 



n satisfying \r\ < s — 1 or \r\ — s, n > k + 1. Pick any particular r satisfying 
we first consider the case of irrelevant vertices. For 2k + z > 5, we may use the same argument 
leading to Eq. Q3.20Q in the proof of Theorem 1 and find that 

pM 



d z v 



A 

r,2k 



(2A,n) 



< 



< 



< 



d z V 



M 
r,2k 



+ 



M 4~2k- 



(2A l( u) 

'Plog(^) 
V m J 



dX 



d x d z V, 



r,2k 



+ 



M 



dX 



(A + mf- 2k - 
W 2 



■ M 

-Plog ( - 



(A + m) 



-2k- 



M 2 



-Plog 



m J 



where, in the second line, we used the induction hypothesis. If 2k + z = 
property that dpV^ 2k (0) = for odd dimensional vertices and therefore 

2n-l 

dt Pi^ d 7 V r,2k(l^ • • • > <l2n-l)\ qi =tp, 



(2A lA t) 



(4.12) 

5, we may use the 
(4.13) 



to obtain the bound (O). The case with relevant vertices 
As in the proof of Theorem |l|, we use the boundary conditions (|4.11| ) at A 

r A 



(2A,M) 

i.e., with 2k + z < 4) is quite trivial. 



and obtain 



\d w V 



•(/)A 



dX 



d x d z V 



A 

r,2k 



(2\,p) 



(4.14) 



Then Taylor's formula (|3.22|) , with C *2 n replaced by T>^ 2n and also using the argument imme- 
diately after Eq. ( |3.22| ), yields the desired bound. | 

Therefore, the difference between amputated connected Green functions of the <f)-ip theory 

and those of the <fi theory is bounded by (j^ m A ~ 2n ~ z Y\og for \pi\ < fi, which is the 
standard statement of the decoupling theorem. Also note that the difference in the vertex 
functions is no longer small if A becomes comparable to the heavy particle mass M, which 
implies that low-energy effective theory is not useful above the heavy particle mass scale. 

In fact, as we discussed in Sec. 1, the decoupling theorem may look more or less obvious 
if one takes the view that only the light field is subject to direct observation. One may then 
imagine a pure 0-field theory with the Lagrangian given from the Lagrangian of the (f)-ip theory 
after integrating out the heavy field, and identify the cutoff Ao with the heavy particle mass 
scale M. What we have shown above is then that (i) integration of the heavy field and the 
"identification" of the cutoff A with the heavy particle mass scale can be achieved by our 
"mass-dependent scheme", Eq. ( |3.5| ), and (ii) the resulting nonlocal irrelevant terms are indeed 
well-behaved so that decoupling theorem may follow. 
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4.2 Factorization of virtual heavy particle effects 

We now examine higher order terms in 1/M 2 , which originate from virtual heavy-particle 
interactions. Let us suppose that we want to describe low-energy physics in the full theory 
accurately up to order 1/M 2N °. (No is some positive integer.) We will show below that we can 
then factorize all heavy-particle effects to the given order by making use of Eq. ( [2.1 0|) (with 
k = 1) and by appropriately choosing the operators L Ao;N , N = 1,2, ... , No. 
Here, L Ao;7V 's may be assumed to have the form 

tAo-n _ [j4 /local, Euclidean invariant, even polynomials \ _ 1 (aik} 

-J dx { in and its derivatives, of dim.< 4 + 2N J' ~ ' ' ' ' ' ' ' [ b) 

The coefficient of the polynomials will be chosen later so that L A °' N carries information appro- 
priate to O (1/M 2 N ) -effects from the full theory. Then the sum of L Ao;Ar, s, viz., 



1 1 
■A ;iVo = l a + J2 J-L A °-' N (= /> ; aw + ——^L Ao ' N °^ , (4.16) 



N=l 



will reproduce the predictions of the full theory at low energy with an accuracy of order l/M 2No . 
The quantity L A ' N , which satisfies the flow equation ( |2.10|) , may be expanded in powers of 0: 



uvj ^ f 2n-l 74 

L A]N [4>] = II Wr^) ■ ~<MP*n)*$Z(Pl, ■ ■ " >P2n-l) • (4-17) 



r=0 n=l j: 



We may also define C r ^nS through the same kind of relation. Then the vertex function C r ' 2n is 
C°° and satisfies the properties similar to those of C^n'i m particular, there is a finite number 
n (r, N) such that £^2n = if n > n (r, N). C r '^s obey the flow equation 

dA^ZiPl, ■ ■ ■ ,P2n-l) 

/2n + 2\ f d 4 q A AN 

= ~y 2 ) J j2^y dAP ^ C rkn+l)^ -q>Pl> ■ ■ ■ >P2n-l) 

+2E e u'hPiiiP) ( N r) [£%i{pu ■ ■ ., P 2i-i)ct;',y(-p, P 2i, . . . ,p 2 n-i)] mm . 

1=0 r '+r"=r ^ ' syrnm 

l + l'=n + l 

= rfiniPl, ■ ■ -lP2n-x) , (4-18) 

where P = —J2f=iPj- To complete the definition of L A °' N , we must also state the boundary 
conditions for C^'^. At A = A , Eq. ( |4.15[ ) implies that 

dp£r$n = 0, for 2n+\w\>A + 2N. (4.19) 
For 2n + \ w\ < 4 + 2N, we impose the following conditions, recursively in N: 

^ST(O) = iV!3;(4gf°(0) - £^n' N -\0)) , (4.20) 
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where £ is defined in Eq. (pT9|) . (Consequently, <9™£ r ^ n ' (0) = for the case 2ra+|w| < A + 2N.) 
It remains to show that the boundary conditions ( |4.19| ) and ( |4.20|) define L Ao;Ar 's uniquely. For 
this, we use Eq. (|4.18| ) together with the conditions (|4.19|) and ( f4.20|) to rewrite 



fA ;N _ 4+2N-2n r-A=0;N 
*-r,2n ' *-r,2n 



Ao 



d\T 4+2N - 2n ^f n , 



(4.21) 



where the symbol r l stands for the Taylor operator defined by 



T l f(pi,...,Pn) 



(for I > 0) 



(4.22) 



t=o 



(and r l f = for I < 0). Then using induction on N, r, and n one may easily confirm that 
L Ao;iV 's are defined uniquely. 

The theory being completely specified, we will now prove the factorization. To begin with, 
we define T> A '^, where r is a multiindex r = (r 1; r 2 , r 3 ) from now on, as 



T^JV _ Af)A 
u r,2n — *- 



-A:N 



~A;N 
— C 
r,2(n,0) ^ r,2n 



(= vXr 1 - c'X/n< 



where C'^ n (and £Q^, correspondingly) is defined by the expansion, = T,r 2 ,r 3 9?9?£'(rur 2 ,r 3 ),2n- 

The flow equation satisfied by D A ^ n wu ^ then be the difference between those given in Eqs. ( p.lOj ) 
and ( |4.18| ). After some algebra, we find that the equation can be written in the form 



(4.23) 



,f2 n rz r/A;N 



o t-.A;V 



2n + 2 
2 



d A q 



(27T) 



<9a P a V A]N + <9a P a f if)A 

<J A^ m ^ ri 2{n+l) + U ^ r M L 'r,{2n,2) 



(4.24) 



+2 £ H^P, 



7" -f-7" =r 
I + !' = n + l 



iV 



El -/A;/ rp.^N-I or /A ^A;N -J-.A ^p.A;N 
T , *~ r',2l U r",2V ^ r' ,21 U r" ,2V U r' ,21 U r" ,2V 



.1=1 



n 



symm 



where, for simplicity, we have omitted the arguments of vertex functions which are the same as 
those in Eq. (|3.10|) or (|4.18 ). The boundary conditions satisfied by 2? r 2n are very simple: 



£nwr\Ao;N |-> 

u p u r,2n — u • 



at A = A : 
at A = : d™V A =°' N (0) = 



for 2n + \w\ > 4 + 2N , 
for 2n+ \w\ < A + 2N. 



(4.25) 



Note that if it were not for the inhomogeneous term containing m Eq- (|4.24 ) , our boundary 
condition Q4.25] ) would imply that T> A '^ should vanish identically. In other words, a nonzero 
result for the difference between the vertex functions of the (p-ip theory and those of the 
theory is due to the inhomogeneous term which is in fact nonzero only for large A, i.e., for 
A > M. As was done in previous sections, we can obtain an estimate for T> 



A;N_ 
r,2n • 



d A d z V A f n 



(2A )A ») 



< const |(A + m) 



d z V 



A;N 

r,2(n+l) 



(2A,At) 



+ #(A-M)A 



«*/•(/) A 
u '~r,(2ri,2) 



(2A,/i,M) 
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+ E 



r -\-r =r 
l + l' = n + l 
z l+ z 2 + z 3= z 



(A + m 



1 / N , 
™) 3+Zl \to 11 



A;I 
2/ 



Z3 r /A;N-I 



(2A,At) 



,2/' 



+ 



(2A,M) 



° u r''2V 



(2A,n) 



(2A, M ) 

. (4.26) 



Now we state the factorization 

leorem 4 The difference betwe> 
of the (p theory, satisfies the bound 



A:N 



Theorem 4 The difference between two vertex functions, Cylti o) °f the cj)-ip theory and C' r 2n 



d z {C 



(/)A 
r,2(n,0) 



~A;A\ 
~ *~ r,2n/ 



< < 



A+m ,2N + 2 

j (A + m) 



A-2n-z 



(2A )A t) 



M 

m) 



Plog 



m) 



< A < M 



A 



4-2n-z 



Plog 



mj 



M < A < A . 

(4.27) 

Proof. Here we use induction on A", r and n. The AT = case has been proved already in 
Theorem || Therefore, supposing that Eq. ( 4.27 ) is true for N < N' — 1 with some N' > 1, we 
may demonstrate the validity of Eq. ( |4.27| ) for A" = N'. Then the rest of the steps are actually 
very similar to those used in proving Theorem |l] or Theorem |3|. If r = or if n > n (\r\, N') 

for r with |r| > 0, we know that £^2(no) anc ^ ^'r,2n (and hence V^'^ also) vanish identically 
and Eq. (|4.27| ) holds trivially for this case. So we prove the validity of Eq. ( |4.27| ) for d^V^k 
with |r| = s under the assumption that it has been proved for |r| < s — 1, and for |r| = s and 
n > k + 1. Under this induction hypothesis, Eq. ( |4.26| ) is now simplified as 



f (A + m) 5+2iV '- 2fc -* , fM\ 
-Plog — 



d K d z V 



A; AT' 
r,2fc 



(2A,/i) 



< ^ 



M 27V '+ 2 

^3+2^-2^-^ 



\m J 



M 2N> 



-Plog 



m/ 



< A < M 



M < A < A , 



(4.28) 



where we used Eq. (^4.23| ) in estimating d Z3 C'^'^. 1 
downward from A = A as before and conclude that 

rAo 



(2A,aO 



. Now if 2k + z > 4 + 2 AT', we go 



d z V, 



A;N 
,2k 



(2A lA i) 



< 



< < 



dX 



d,d z V 



\;N' 
;2k 



(2A,m) 



f (A + m) 6+2N '~ 2h ~ z (M\ 
M 2N '+ 2 ° g \m) 

pA+2N'-2k-z 

M 2 ^' P1 ° g U ' 



< A < M 



M < A < A . 



(4.29) 



(If 2A; + z = 5 + 2N', the upper inequality of Eq. (jOg) is obtained using d^V^' (0) = 0, as 
was done in the proof of Theorem [|) For 2 A; + z < 4 + 2 AT', we have 



(0)| < / d\ 



d.d z vt\ 



\;N' 
2k 



(2A,m) 



(4.30) 
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and so Eq. (|4.27| ) follows with the left hand side replaced by \d^T>^^ (0)|. Then, again, Taylor's 

formula (|3.22|) (now applied to T)^'^. ) relates the corresponding vertices with those involving 
twice more differentiation by momenta. Therefore, to complete the proof, induction on z should 
also be used (for fixed k). If z > 4 + 2N' — 2k, Eq. ( |4.27| ) holds as we have just shown above. 
Now assume that Eq. (|4.27 ) is true for all z > z for some z > 0. Then, finally, Taylor's 



formula and Eq. ( (4.30| ) guarantee the validity of Eq. ( f4.27| ) for any w with \w\ = zq. I 



If we choose A = 0, in Eq. ( |4.27| ), it yields the bound for the difference between connected 
amputated Green functions of the (p-ip theory and those of the <fi theory with insertions of 
appropriate operators defined above. Explicitly, for \pi\ < \i, 

\d;(G { r % >0) ( Pl , . . . ,p 2n _ 1 ) - Gf£°(Pi, ■ ■ ■ ,P2n-i))\ < {^) 2N0+2 m '~ 2n ' Z P1 °S (^) > ( 431 ) 

where G^n — £-2n° ,N ° = ^n°=o G^ /N\. Thus one may identify the effective theory which 
is accurate up to order l/M 2N ° as that defined by L A °' No = L A ° + T,n°=i L Ao ' n /N\ and the 



flow equation ( 2.10|) . We will discuss more on this effective action in Sec. 4.3. If A is larger 
than M, Theorem (| says that the vertex function dpC^2n is larger than dpCf,'2n~° by a factor 
(A/M) 2JV . Therefore, in this scale, increasing (or adding more irrelevant terms, equivalently) 
makes the vertex functions "unnaturally" large and the theory becomes "less effective" . This 
is because we have imposed unnatural renormalization conditions ( |4.2U| ) on irrelevant vertices. 
If we chose natural values as the bare irrelevant couplings, they would give only O((m/A ) 2N ) 
contributions to Green functions; but, in Eq. (|4.20|) , we imposed the values of order (m/M) 2N 
to CI ' 2n as the renormalization conditions, which are natural only if the cutoff is around M. 
This affirms that the effective theory is useful only below the heavy particle mass scale M. 
The convergence bound as A — > oo shows a similar behavior. For A < M, the difference 

converges faster by a factor i^jf 1 ) than the A^ = case, while the convergence becomes 
worse if A > M. We state the result here omitting proof: 



-A;A\ 
r,2n / 



s {—) {—) (A + m) Pl0g UrJ • (432) 



4.3 Irrelevant operators in the effective Lagrangian 

So far we have shown that, at low energy, connected amputated Green functions of the 
full theory are reproduced by considering those of the <fi theory plus those with insertion of 
operators On, N = 1,---,N , defined in Eq. (|2.8| ). Such operator insertions have a simple 
interpretation if one uses the normal product notation |15|, || described in the Appendix. Let 
us begin with Ao = 1. Comparing Eq. ( |4.20| ) with Eq. ( |A.3j ), we identify L Ao;1 as 



> !l = E b% {w} [M 2n , {w} ] , (4.33) 

2n+\w\=6 



18 



where = jr^^p(^2(no) ~ ^2n)(0)> -^{w} is a combinatorial factor defined by Eq. ( |A.3| ), 

and [M 2ni { w }] is the normal product of a local vertex 



M 2n ,{ w } 



= Jd^d^^-'-d^-'^ix) (4.34) 
(see Eq. ( |A.4| )). In terms of the local Euclidean-invariant vertices of dimension six, the above 



equation can be cast into a more attractive form given in Ref. |13 



j>i = a ? ! d 4 x [0(<9 2 ) 2 0(x)] + af ( d 4 x [0 3 <9 2 0(x)] + af [ d*x [0 6 (x)] 



= E«fN 6) ], (4-35) 



with 



a 



(6) _ 1 (^2 (n (f)c 



(9 2 ) 2 (g;5-g c 2 )(o) 



8-4! 

T (4,0) 



4 6) = -^(G ( &-Gt)(0) 
4 6) = (G$J-G§)(0). (4.36) 

The l/M 2 -order effects are then described by the insertion of the operator 0\ = J2i a f^ [^f^]'- 

3 ... / m &-2n^ fM , - 

m. 



G%&) = G ^ + E o^G£{[^)) + O [ ^Plog ( ^ ) ] , (4.37) 



i=i 

where G ? 2n([^' 6 ' ) ]) denotes the Green function with \ inserted. For general ATq, L A °' N (N 
1, . . . , N ) can be identified as 



l a ; N= ^ b«% w} [M 2n , {w} ) + L% N , (4-38) 

2n+\w\=4+2N 

where 

b 2n,{w} ~ Tf— 

- /v {w} 2n+\w\=4+2N 



W:Y E d^G^-G^m, (4.39) 



and Lqj? 1 denote counterterms which are needed to cancel the new divergences due to the 
multiple insertion of L Ao;/ 's, / = 1, . . ., N — 1. (See the discussion of the last paragraph in 
the Appendix for explicit forms of L^ N .) Now let [f2- 4+2Ar ' ) ]'s form a complete set of mutu- 
ally independent Euclidean-invariant local vertices of dimension (4 + 2N) and let a\ 4+2N ^ 's be 
appropriate expansion coefficients as in Eq. ( EL35D , so that we may write 

L A °' N = £ at 2N) [ff +2A °] + L%f . (4.40) 
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Then the effective action L A °' No may be written as 

No J4+27V) iVo -r 

~ L A ;No = LHEEV^ 1+^1^ 

N=l i JV - N=2 JV ' 

= Z a ;^o-i + v^ [o( 4 + 2 ^)] + -Ll^ , (4.41) 

where a[ 4+2N ^s are appropriate linear combinations of b^fi w y8 (all of which are 0(1/M 2N )). 

Thus the l/M 2Ar -order information of the full theory is factorized into the coefficient s of 

local vertices of dimension (4 + 2iV). Given this bare Lagrangian, the flow equation guarantees 
that all Green functions are finite and accurate up to order 1/M 2N . 

The interpretation of this result is clear from the viewpoint of the renormalization group 
flow in the infinite dimensional space of possible Lagrangians. As we reduce the cutoff, the bare 
Lagrangian of the full theory, which is specified by the boundary conditions ( f3. 11 ) and ( 3.12j) , 



flows down to a submanifold parametrized by relevant couplings only. It has both heavy and 
light operators. Now the flow may be projected onto the submanifold of operators consisting 
of the light field only. Then finding a local low-energy effective field theory of light particles 
is equivalent to finding a renormalization group flow with a few (relevant or irrelevant) local 
operators which can best approximate the projected flow of the full theory. To the zeroth order 
in 1/M 2 , it is done with relevant terms only by choosing the same renormalization conditions 
as those of the full theory. The deviation from the full-theory flow is corrected by reading 
off values of the remaining irrelevant coordinates at the A = point. At first, components of 
dimension six operators may be read, which tell us the effects of order 1/M 2 ; the flow of the full 
theory is now approximated to the order 1/M 2 at A < M. If one wants to increase the accuracy, 
it is necessary to read more and more irrelevant coordinates (at A = 0) of the projected flow of 
the full theory and modify the effective-theory flow with the corresponding flow equation given 
in Eq. ( |2.10| ). Our equation fl2.10| ) automatically takes care of possible divergences due to the 



unnaturally large irrelevant components (see the discussion of Sec. 1) in such a way that the 
bare Lagrangian may contain only a finite number of irrelevant terms. 



5 CONCLUSIONS 

In the context of a simple scalar field theory, we have demonstrated that, at low energy, 
virtual heavy particle effects on light-particle Green functions are completely factorized via 
effective local vertices to any desired order. For this, we have used the powerful flow equation 
approach which can be regarded as a differential version of Wilson's renormalization group 
transformation. Remarkably, the method which has been used previously to prove renormaliz- 
ability is found to have an immediate generalization to this problem. In applying this method, 
we have seen that irrelevant terms in the Wilsonian effective Lagrangian (which has a finite 
UV cutoff A = M, a characteristic scale representing heavy-particle thresholds) are replace- 
able by the corresponding higher dimensional composite operators plus counterterms for their 
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products in the conventional approach (where UV cutoff is supposed to go to infinity). The 
latter can be dealt with with the help of the normal product concept. Once this fact is no- 
ticed, factorization is straightforward: all l/M 2JV -order effects can be isolated in terms of local 
vertices involving operators of dimension (4 + 2N). Thereby we arrive at a local effective field 
theory which describes low-energy light-particle physics accurately up to any desired order in 
1/M 2 with appropriate calculation rules for irrelevant (nonrenormalizable) vertices. Since the 
arguments here are essentially dimensional, it should not be difficult to generalize them to 
different field theories such as gauge theories or theories with spontaneous symmetry breaking. 
As for these models, we note that the flow equation approach has already been used to discuss 
renormalizability and related problems P, ^, £|. 

In this paper we limited our attention to vertex functions at momentum range \i = 0(m). As 
one increases \i to sufficiently high energy scale, the (/x/m)-dependence we neglected will become 
important. Keeping such /^-dependence and establishing bounds more carefully, one may study 
the high-momentum behavior of Green functions (even in the exceptional momentum region) 
within the effective field theory context. In this regard, we note that the infrared behavior 
of massless scalar theory has been analyzed by using flow equations recently [^, ||. Another 
related point is the following. Effective theory is supposed to be useful mainly at low energy. As 
we increase energy, we need in fact more and more irrelevant terms to maintain the same level 
of accuracy. Conversely, with a few irrelevant terms added to the Lagrangian, we may extend 
the energy range where it remains effective. Then it seems worth investigating quantitatively 
to what energy scale the given effective theory may be used to reproduce the full-theory results 
with reasonable accuracy. 



NOTE ADDED 

Completing the main parts of this paper, we became aware of the very recent paper by 



Girardello and Zaffaroni [20] which has a partial overlap with the present paper. They also 
discussed decoupling and factorization to order 1/M 2 using flow equations. However, we believe 
that our treatment is physically more transparent and systematic. For example, they did not 
discuss the normal product description nor higher order generalization; this does not look trivial 
in their method for, in their approach to these issues, nonnegligible parts of heavy-field modes 
might remain unintegrated even at scales much less than the heavy field mass M. 
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APPENDIX 



Here we briefly discuss the renormalization of composite operators and their products |4], [|, 
within the 4 -theory defined by L A ° in Eq. (|4.1|) . We restrict our discussions to those needed 
in the main part of this paper. The operator L Ao;Ar, s of Sec. 2 are assumed to be of local, even 
polynomials in </> and its derivatives of dim.< 2D^, 



L A ° 



D N 

N = H E 

n=l 2n+\w\<2D N 



^2n,{w}^2n,{w} , 



(A.1) 



where b* n{w} 



is a formal power series in g 1; 6^ r , 



1 '' A which is to be determined 



later uniquely when we impose renormalization conditions. [M 2n ,{w} is defined in Eq. ( f4.34j) .1 
Throughout this Appendix, the dimension is restricted to satisfy the condition (2D^ — 4) > 
Ei=i N i( 2D i - 4 )' where A = D N \ N . =S .., N = (N u ...,N k ). We follow the now familiar 



procedure. £^ 2n' s are defined by Eq. (|4.17|) and have similar properties except for the Euclidean 



invarianceQ. They satisfy the flow equation ( |4.18| ). Boundary conditions are given similarly: 



A = A : 
A = : 



d w p L 



A ;N 
r,2n 



h R;N K r 

U r;2n,{w} J ^ {w} 



for 2n + \w\ > 2D N , (A.2) 
for 2n+ \w\ < 2D N , (A.3) 



where M{ w ] is a combinatorial factor defined by d™' [(ipi] 



Wl 



(ip2n-l) W2n - 



i symm 



j ,R;N 
and K;2n.{ W } 
U R;N 




p=0 



S 



{■w}{w'}N{ w} ) 



s are finite and A -independent constants. [Here it is to be noted that b^ n r tt i and 



2n {w} are in g enera l n °t the same; this is the usual operator mixing 0, |].] With these one can 
easily prove the perturbative renormalizability of Green functions with the insertion of operator 
On (which is defined in Eq. ( |2.8| )); the proof is just a simpler version of that in Theorem 4. 

We now discuss Zimmermann's normal product following Ref. ||[]. We begin with a 
monomial for the case of k = N = 1 with D% = D. Normal product [Mm' ,{w'}]2D of dimension 
2D (> 2n'+\w'\) is defined by requiring that the vertices with one insertion of 0\ = [M2 n >,{w'}]2D, 
which are denoted as C r 2 n {[M2n' ,{w'}]2d) i obey the renormalization condition (|A.3|) with 



,R:1 



u r;2n,{w} ~ $nn'${w}{w'}$rO ■ (A .4) 

Since the flow equation is linear, for general boundary condition (|A.3| ) we may write L Ao;1 as 



L A ; 



E 



■ R:l 



2n+{w}<2D 



hn,{w}l M 2n,{ 



\2D 



E 

2n+{w}<2D 



b 2n,{w} M 2n,{w} 



(A.5) 



2D 



5 Since the flow equation respects 0(4) Euclidean symmetry, one may well restrict the consideration to 
Euclidean- invariant operators and then C r ^ n will also be Euclidean invariant. 

6 Here we are discussing only space-integrated operators, or vertices. Generalization to local operators, e.g. 
[M 2n ,{w}(x)} 1 are straightforward and is described in Ref. || to the order k = 2 and \N\ = 2. 
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When the dimension of the operator is equal to 2D, we simply write [• ■ -] 2 d = [•••]; f° r example, 
[M 2 n t { w }]2D = [^2n,{w}] if 2D = 2n + \w\. Therefore, if we choose hf^^y = for 2n+ \w\ ^ 2D, 
we may omit the subscript D in Eq. ( |A.5|) . This is the case of Sec. 4 (see Eq. Q4.33j )). 



The normal product of several operators (assuming k > 1) are defined inductively as follows. 
Let us write N = L A °'> N = [B^. for N = (N u . . . , N k ) with Nj = <J y (i.e., \N\ = 1). 
Suppose that [[-B il ] 2 D il • • • [Bi k _j2D ifc _j2D's have been defined, where (2D— 4) > Ef=i 1 (2D ii -4), 
1 < ij < Then we define the normal product [[-Bi] 2 Di • ■ • [Bk\2D k ]2D N as 0(1,...,!) which flows 
down through Eq. (]2.10|) with the simplest boundary condition 



^/^5? (W) (0) = b^^Af {w} = , 2n + | W | < 2D,v . (A.6) 

If (2Dtv — 4) = X)i(2Dj — 4), we write [• ■ -]2D N = [•••] as above. Let us illustrate this definition 
for k = 2 and k = 3 cases. If k = 2, [[-Bi^x^ [D 2 ]2d 2 ]2D (1 1} is defined as 

[[Bi] 2Dl [B 2 ) 2D2 ] 2D(1A) = (1>1) = [D x ] 20l [D 2 ] 2D2 + L Ao ^ , (A.7) 



where the flow of L^ '^ 1 ' 1 ' is governed by Eq. (|2.10 ) together with the boundary condition ( A.6 ). 



Therefore, D^ 1 ' 1 ) is interpreted as the counterterm needed to cancel divergences appearing in 
Green functions with the insertion of the product operator [Bi] 2 d 1 [B 2 ] 2 d 2 . Similarly, if k = 3, 
from the expression for 0(1,1,1) given in Eq. ( |2.9|) we have 



[[Di] 2Di [D 2 ] 2jD2 [D 3 ] 2D 3] 2D(1i1i:l) = 0(1,1,1) (A. 8) 

= [Di] 2i?1 [B 2 \ 2D2 [B 3 ] 2D3 + [B^L^ ' 1 ^ 

+ [B 2 ] 2D2 L A °^ + [B 3 ] 2D3 L Ao ^ + L A °^ , 

where D^oa 1 ' 1 ' ) etc. are given by Eq. (|A.7|) , and additional divergences due to multiple insertions 
are cancelled by L Ao;< - 1,1 ' 1 ' ) which obeys the flow equation ( |2.1(J| ) with iV = (1,1,1), with the 
boundary condition (|A.6|) . In this way, [[Di] 2 d 1 • • • [Bk] 2 D h ]2D N is defined for general k, i.e., 
^a ;(i,...,i) serves as "counterterm" for new divergences appearing in Green functions due 
to the new products of lower order operators and [[Di] 2Dl • • • [B k ] 2Dk ] 2DN is interpreted as the 
resulting "subtracted" finite part. Since this definition satisfies multilinearity, which can be 
proved by use of induction on k while taking into account Eqs. (2J3) and (|2.10| ), and boundary 



conditions, we can freely put the coefficients multiplying operators inside or outside the brackets. 

For the general boundary conditions ( |A.3| ) where b 2 ^i w \ need not be zero, we recall the 
property of the flow equation ( |2.10| ) mentioned at the end of Sec. 2. The general solution of the 



flow equation is represented as a sum of a particular solution and a solution of the homogeneous 
equation which is identical to the flow equation with \N\ = 1, i.e., 

L a -' n = L a -> 1 ([B n } 2Dn )+L^. (A.9) 

Here we may absorb all the nonzero boundary values into L Ao;1 ([Djv] 2 Div) choosing 

[B n ] 2Dn = E b* N {w} [M 2n>{w} ] 2DN , (A.10) 

2n+\w\<2D N 
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and then satisfies the simplest boundary condition ( |A.6| ). (The subscript CT is attached 



to indicate its counterterm nature as discussed above.) In other words, &^fLa 7^ corresponds 
to adding new composite operator [B n ] 2 d n to L A °' N besides the counterterms needed just to 
cancel divergences associated with the multiple insertion of operators already included. For 
example, if k = 2, we may extend Eq. ( |A.7|) to the general case, 

L Ao;(i,i) = L W) + [ B{hl) } 2D(11) (A.ll) 

= [[Bl}2D 1 [B2\2D 2 }2D {1X) - [Bl]2D 1 [ B 2\2D 2 + b 2nf{w}[ M 2n,{w}]2D (hl) ■ 

2n+|to|<2_D {11) 

This is the case we encountered in Sec. 4. 
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Figure caption 



Fig. 1. Graphical representation of the flow equation ( 3.10|) . Thin (thick) lines represents 
light (heavy) field, and numerical factors in Eq. (|3.10| ) are ignored here. 
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